Possible solution of the Coriolis attenuation problem 



Pavlos Protopapas * and Abraham Klein t 
Department of Physics, University of Pennsylvania, Philadelphia, PA 19104-6396 

(February 9, 2008) 



Abstract 

The most consistently useful simple model for the study of odd deformed 
nuclei, the particle-rotor model (strong coupling limit of the core-particle 
coupling model) has nevertheless been beset by a long-standing problem: It 
is necessary in many cases to introduce an ad hoc parameter that reduces 
the size of the Coriolis interaction coupling the collective and single-particle 
motions. Of the numerous suggestions put forward for the origin of this sup- 
plementary interaction, none of those actually tested by calculations has been 
accepted as the solution of the problem. In this paper we seek a solution of 
the difficulty within the framework of a general formalism that starts from 
the spherical shell model and is capable of treating an arbitrary linear combi- 
nation of multipole and pairing forces. With the restriction of the interaction 
to the familiar sum of a quadrupole multipole force and a monopole pairing 
force, we have previously studied a semi-microscopic version of the formalism 
whose framework is nevertheless more comprehensive than any previously ap- 
plied to the problem. We obtained solutions for low-lying bands of several 
strongly deformed odd rare earth nuclei and found good agreement with ex- 
periment, except for an exaggerated staggering of levels for K = ^ bands, 
which can be understood as a manifestation of the Coriolis attenuation prob- 
lem. We argue that within the formalism utilized, the only way to improve 
the physics is to add interactions to the model Hamiltonian. We verify that 
by adding a magnetic dipole interaction of essentially fixed strength, we can 
fit the K = i bands without destroying the agreement with other bands. In 
addition we show that our solution also fits 163 Er, a classic test case of Coriolis 
attenuation that we had not previously studied. 
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I. INTRODUCTION 

In this paper we propose a solution to a long-standing conundrum in the theoretical 
description of the properties of low-lying bands of deformed odd nuclei within the framework 
of the particle-rotor model Generally speaking, the Hamiltonian of this model has the 
form 

H = H TotoI + -^particle- (1) 

Here the first term is a phenomenological operator describing the included bands of a given 
even (core) nucleus, assumed to be a rotor. In the simplest case, it describes only the ground 
band of an axially symmetric system. The second term contains the single-particle energies 
of and the residual interactions among the particles that move outside the rotor and are 
assumed to be kinematically independent of the constituents of this core. When we wish to 
describe odd nuclei that are adjacent to the given core, the simplest assumption is that it 
suffices to consider a single extra particle, or if we include pairing interactions, quasiparticle, 
outside the core, described in the intrinsic system of the rotor by a set of axially symmetric 
solutions of the Nilsson-BCS formalism. 

An important talking point in favor of this extreme model is that it contains, in so 
far as its predictions of energy levels is concerned, essentially no free parameters, since 
the quasiparticle properties can be determined by fitting bandheads alone, and the core 
Hamiltonian depends to a good approximation only on the observed moment of inertia of 
its ground-state band. When the angular momentum of the core is replaced, according to 
angular momentum conservation, by the difference between the total angular momentum J 
and the particle angular momentum j, the interaction between particle and core becomes 
(omitting BCS occupation factors) 

LA 

21 ' 

where X is the moment of inertia of the core. 

For many applications, the model as described above does not predict the observed band 
structure |2||| unless one multiplies the Coriolis interaction by a factor £ = .6 — .8, thus 
giving rise to the name Coriolis Attenuation. This means that the model as described above 
is too rudimentary, and that there is missing physics. There have been numerous suggestions 
concerning the source of this missing physics. We shall not review these suggestions in the 
form that they are described in the literature for at least two reasons. The first is 
that all those that have actually been tested by calculations fall short in some respect. The 
second is that they are all stated within the framework of some extended version of the 
particle-rotor model. But this model has among its flaws the assumption of the kinematic 
independence of the particle and rotor variables. This certainly violates the Pauli principle, 
the more so if one wishes to allow the particle to have access to several major shells (as in 
the calculations we report in the present work). 

Instead we shall approach the subject from the point of view of a theoretical framework 
that is formally exact and contains the particle-rotor model as a well-defined limit. We 
claim that for a sufficiently simple Hamiltonian we can, within this framework, compute 
accurate solutions for the low lying bands of well-deformed odd nuclei. Therefore with high 
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probability a failure of the theory is to be ascribed to a missing term in the Hamiltonian. 
For example, one of the classes of proposed explanation within the particle-rotor framework, 
that we need more than one particle outside the core, is ruled out in our case by the fact 
that the exact formalism never involves more than one particle outside a band of the core, 
the point being that the inclusion of core excitation is taken care of by the inclusion of 
excited bands of the core. Another explanation of the standard list is that the model in its 
simplest form omits or treats incorrectly the so-called recoil term that arises from the same 
transformation that yields the Coriolis coupling; the distinction between the two terms is 
a consequence of the use of an intrinsic frame in the implementation of the model. But 
our calculations are carried in the so-called laboratory system, in a way that automatically 
includes the effects of the recoil term. The last important class of suggestions and the 
only one that is cogent from the point of view of the theory we utilize is that of a missing 
interaction. Here, the structure of our equations suggests an essentially unique choice, at 
least as far as multipolarity is concerned. Further discussion of this point, the major one 
made in this paper, is continued below, after some remarks about the formalism. 

The theory of collective motion for odd nuclei that we shall apply is one that we have 
recently resuscitated and extended, and that we have rechristened the Kerman-Klein-Donau- 
Frauendorf (KKDF) model PHSj. This model was introduced by Donau and Frauendorf 
P-I^l^f, whose work was in turn stimulated by an application [13| of the theory of collective 
motion developed by Kerman and Klein |14|-|17| • This theory starts from the spherical shell 
model with residual two-body interactions. A formally exact set of equations is derived that 
resembles the equations of a particle-core coupling model, but requires even after a suitable 
restriction on the size of the configuration space, the solution of a formidable non-linear 
problem. In the KKDF model the problem is linearized by choosing the matrix elements of 
multipole and pairing operators, that arise from the interaction, either phenomenologically 
or from experiment. If this can be done for a sufficiently extensive model space, we can be 
reasonably confident that the corresponding equations of the non-linear theory have been 
solved accurately (although we have no proof of this statement). The KKDF program can be 
carried out at the present time only if the interaction chosen is sufficiently simple that in fact 
the "core" matrix elements that enter can be determined at least in part from experiment. 

In previous work we have used a conventional Hamiltonian in which the interactions 
were confined to monopole pairing and quadrupole multipole interactions. Among the ap- 
plications made, the simplest have been to the low lying bands of several well-deformed odd 
rare earth nuclei. To a good approximation the only free parameter was the strength of the 
quadrupole interaction. With the important exception to which this paper is devoted, we 
found good agreement with the energies and electromagnetic moments (where measured) 
of all the lowest lying bands, provided the calculations were carried out with three major 
shells for the single-particle space and with the inclusion of excited bands of the core (in- 
cluding but not confined to beta and gamma bands). The inclusion of the excited core 
bands was necessary to account for odd bands near the top of the energy region considered, 
approximately 1 Mev. 

The only discrepancy found was that our calculations predicted, uniformly, an exagger- 
ated staggering of levels for K = | bands compared to experiment. In a particle-rotor model, 
it is easy to see from standard formulas U that this discrepancy can be erased by introduc- 
ing the Coriolis attenuation factor. As we have argued above, the only fundamental way 
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that is open to us to improve our calculations is to add an interaction to the Hamiltonian. 
Such an interaction should effectively contribute to our equations a Coriolis-like coupling 
that is dynamically independent of the standard one (the latter hidden in our formalism 
because we calculate in the laboratory system). Once we have reached this point in our 
thinking, it is almost obvious that the choice should be a magnetic dipole interaction, also 
satisfying our requirement that the core matrix elements can be measured. In fact, from the 
well-known properties of vector operators it follows that contributions diagonal in both core 
angular momentum and particle angular momentum of any dipole interaction are exactly 
equivalent to a Coriolis interaction. However, the magnetic dipole interaction can also have 
off-diagonal elements and therefore the equivalence is only approximate. 

It may be noted that there is even a simpler way to obtain an additional Coriolis coupling 
than the one we have chosen, and that is to add to the Hamiltonian a term proportional 
to the square of the angular momentum. We shall verify that such an interaction can also 
be used to fit the data. Such a Coriolis coupling term was seen to arise, for example, in 
a paper [18[ showing that a version of the particle-rotor model can be derived from the 



Fermion Dynamical Symmetry Model (FDSM). This is not surprising, since the standard 
Hamiltonian of FDSM indeed contains a term proportional to the square of the angular 
momentum operator of that model. This observation was not followed up. 



It is also of interest to remark that in a recent analysis |19| which sought to identify 
the most important multipole components of a shell-model interaction, it was found that 
magnetic dipole interactions were present. The result was represented as a spin-spin inter- 
action and therefore different in detail from either of the possibilities examined in this work. 
Nevertheless this analysis further buttresses our viewpoint. 

Finally we outline the contents of the body of work that follows. In Sec. II, we briefly 
review the basic elements of our theory and incorporate the changes necessary to include 
the new interaction. In Sec. Ill, we specify how the parameters of the magnetic dipole 
interactions were chosen. The results of our calculations are presented in Sec. IV, together 
with a brief description of how the analysis was carried out, where it was deemed necessary 
to supplement discussions in our earlier papers. In Sec. V we compare some results found 
for the magnetic dipole interaction with a corresponding analysis in which we simply add 
to the Hamiltonian a term proportional to the square of the angular momentum. A brief 
discussion of our results is given in Sec. VI. 



II. THEORY 



For the sake of completeness, we start with a brief review of the basic elements of our 
method. We base our considerations on a shell-model Hamiltonian of the form, 

H = J2 h a a f a a a + ^ E E F acdb{L) B[ ML (a, c)B LML (d, b) 

a 4 abcdLM L 

+ 7 E E G aUL) 4m>, b)A LML (c, d). (3) 

abcdLM L 

Here h a (a = (j a , m a ) and a = j a ) are the spherical single particle energies, -Blm is the 
elementary particle-hole multipole operator, 



4 



B LM L {a, b )= E s ? C tf L a « a P = [ a l x ^b\ L Mh , ( 4 ) 



m a m b 



and Al M is the elementary particle-particle multipole operator 



4m>, 6) = E Ci Mi 4 4 = [4 x ajf , (5) 



m a m b 



where C^g* are the Clebsch-Gordan coefficients and s a = (— ma . The coefficients F are 
the particle-hole matrix elements 

F acd6 (L) = J2 W C ^ V afhS , (6) 

m's 

and G the particle-particle matrix elements 

G ahcd (L) = V af3l& . (7) 

m's 

The task is to obtain equations for the states and energies of the odd nucleus assuming 
that the properties of neighboring even nuclei are known. The states of the odd nucleus 
(particle number N) are designated as | Jfiu) where v denotes all quantum numbers beside 
the angular momentum J and its projection fi. The eigenstates and eigenvalues of the neigh- 
boring even nuclei with particle numbers (iV±l) are \IMn (N ± 1)) and Ef^ 1 , respectively, 
where n plays the same role for even nuclei as v does for the odd nuclei. The equations of 
motion (EOM) are obtained by forming commutators between the Hamiltonian and single 
fermion operators, leading to 

[a a , H] = h a a a + j E C^fG acM {L)a\ A LM {c, d) 

4 bdj LM 

+ 7 E E s,C^F acdb (L)a,B LM (d, b), (8) 

4 bdj LM 

together with its hermitian conjugate. 

The matrix elements of these equations provide expressions that determine the single- 
particle coefficients of fractional parentage (cfp), 

Vj^(a; IMn) = (Jfii/\a a \I Mn(A + 1)), (9) 
Uj^(a; IMn) = (J H4|/ Mn{A - 1)). (10) 

In terms of a convenient and physically meaningful set of energy differences and sets of 
multipole fields and pairing fields defined below, we obtain generalized matrix equations of 
the Hartree-Bogoliubov form 

£jvV Jlu ,(a; IMn) = (e + + r^ +1 )) a/Mni7/ , M ,„,lV( 7 ; I'M'n') 

+A aIMn ,-fi'M'n'Uj IJ , l/ (Y, I'M'n'), (11) 
SjuUj^a- IMn) = (-e + - ^ A -^) aIMn ^ rM , n ,Uj^(r, I'M'n') 

Vj^ud; I'M'n'), (12) 
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where 



talMn^I'M'n' 



-Ejy + ~(Eq 



(A+l) 



S a ^Su'SMM'^nn'(h a — Xa), 



CO 



(Ail) 

aINn,-yI'M'n' 
aIMn,"/I'M'n' 



*-aIMn,jI'M'n' 



(Ail) 



E 



(Ail) 



), 



L t W {A+l) F (A-1) 
Sa^II'^MM'^nn 

Z LM L bd 

(I'M'n'(A±l)\B LML (db)\IMn(A±l)) 

o EI ^2 C a^ L G acdb {L) 
Z LM L bd 

(I'M'n\A - l)\A LML (db)\IMn(A + 1)). 



(13) 
(14) 
(15) 
(16) 

(17) 
(18) 



Here 8 is the eigenvalue, the energy of the state of an odd nucleus relative to the average 
ground state energy of the neighboring even ones. The remaining definitions refer to the 
elements of the effective Hamiltonian matrix: e are the single particle energies measured rel- 
ative to the chemical potential, A the chemical potential, the elements of u are the excitation 
energies in the even nuclei, and Y and A are the multipole and pair fields, respectively 

The solutions also require a normalization condition that is derived from matrix elements 
of the fundamental anticommutation relations, 
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aIMn 



]T \\U^(aIM;n)\ 2 + \V^(aIM;n) 



where 



n = £(2j. + i). 



(19) 



(20) 
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Since we shall use only a restricted set of interactions, specializing to the case where 
L — 1,2 multipoles and the L = pairing are included, the multipole and pairing fields 
defined above take the form 



,(Ail) 

alMn^I 1 M'n 1 



1 



*-aIMn,jI'M'n' — 



" Mi bd 

(I'M'n'(A ± l)\B 1Ml (db)\IMn(A ± 1)) 

Z M 2 bd 

(I'M'n'(A ± l)\B 2ML (db)\IMn(A ± 1)) 

9 EZ c ™G acdb (0) 
z bd 

(I'M'n'(A - l)\A 00 (db)\IMn(A + 1)). 



(21) 



(22) 



The further specialization to the magnetic dipole, the mass quadrupole, and the standard 
monopole pairing interactions, is defined by the equations 
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F acdb (l) = -K 1 M ab (l)M dc (l), (23) 
F acdb (2) = -K 2 F ab (2)F dc (2), (24) 
G acdb (0) = -g o G ab (0)G dc {0), (25) 

together with the identification of the sums 

M^ = J2M ac (l)B lMl (a,c), (26) 

ac 

Q, = J2 F ac(2)B2M 2 (a,c), (27) 

ac 

A = J2G ac (0)A 00 (a,c), (28) 

ac 

as the electromagnetic dipole, the mass quadrupole, and the pairing gap operators, respec- 
tively. 

With the definitions given above, the multipole and pairing fields take the form 

L alMn^I'M'n' — 9 A 7 ^07 lvI acK L ) 

Z Mi 



x(I'M'n'(A±l)\M Ml (db)\IMn(A±l)} 

K2 

2 



M 2 

x(I'M'n'(A±l)\Q M2 (db)\IMn(A±l)), (29) 

AaiMn^I'M'n' = ~^G^G ac (0) 

x(I'M'n'(A-l)\A\IMn(A+l)}. (30) 



We apply the Wigner-Eckart theorem to obtain the equations of motion in their final 
form. For the multipole and pairing fields, we use the following definitions of reduced matrix 
elements 



(I'M'n'(A±l)\M Ml \IMn(A±l)) = 1 ^_ C)$ V _ M , {I'n\A±l)\\M\\In{A±l)), (31) 



(I'M'n'(A±l)\Q M AlMn(A±l)) = { ^_ C f 2 ^,-M' (/V(A± l)||Q||/n(A± 1)>, (32) 
(J'MV(A±l)|A|/Mn(A±l)) = (-l) J - Af C^ 7 ,_ M ,(/W(A± l)||A||/n(A± 1)), (33) 
and for the cfp V and [/, we employ the definitions 

Vj, u (a, IMK) = ^== C? M ^ V^a/X), (34) 

C/^(a, IMK) = 1 g/Mj-, Wj„(a//0- (35) 

We thus find the equations of motion 
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£ Jv V Ju (aIK) = (e„ + u IK ) V Ju (aIK) 

+ E r&gU VjM'k') 

d'K' 

+ X! Aa/x, c /^' U Jv {cI'K'), (36) 
^ U Ju (aIK) = (-e a + a;/^) U Jv (aIK) 

+ £ Aaf^^ VMcI'K'), (37) 



where 

rif^ = -y(-l) Jc+/+J {^ c j} (/^ / (A±l)||X||/X(A±l))M ac (l) 

- y(-l) ie+/+J {^ j j} (/'^(A ± 1)| jQ|]7ir(A ± l)>i^ c (2), (38) 



= -|(/V(A± l)||A||/n(A± 1)). (39) 



III. DETERMINATION OF MAGNETIC PARAMETERS 

In this section we describe how to extract the magnetic parameters for the core and for 
the single particle motion. 

A. Phenomenology of the even core 

The matrix elements of the operator M. between states of the even nuclei have to be 
determined. The reduced matrix elements of this operator can be expressed using the 
Wigner-Eckart theorem as 



(I\\M\\I) = ?-^(II\Mo\II) (40) 

From the definition of the intrinsic magnetic moment /i we have 

H = (II\M \II) (41) 

Consequently we have 



(I\\M\\r) = ^ (42) 

(43) 

Because of the limitation on available experimental data for the intrinsic magnetic mo- 
ments of the states of the neighboring even nuclei, we have chosen to represent their values 
by the simplest possible phenomenological model, 
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H = g R (I)I ~ gj, 



(44) 



where the parameter g^ is the effective g factor for the rotational motion. Though this 
quantity can itself depend on the angular momentum, we have found that the use of a 
constant value, gi, provides a reasonable, though not perfect, fit to the data and, at the 
same time, simplifies our study of the odd systems. For the resulting values of g\ see 
Table PL 



B. Matrix elements of the Ml operator for a single particle 

For the matrix elements of the Ml operator between single particle states we quote the 
formula from (111 as 



0. II M || i) = (».U | u Mc ) f (2j ° + ^ + 1) (-l)>-'^ ( _{" /2 I l< 2 



x(l — k) 



-g s -gi(l + 2 



where the first factor is a radial overlap integral, and 



k = (ja + + (j c + l)(-l)(Wl/2)_ 



(45) 



(46) 



Here g s and gi are the spin and orbital gyromagnetic ratios. In our calculations, we have 
adopted the values suggested by the analysis found in EQ|, namely 0.70 of their bare values. 



IV. APPLICATIONS 

We have applied the theory described above to four nuclei, the odd neutron nuclei 157 Gd, 
159 Dy, 163 Er and the odd proton case 157 Tb. All but 163 Er were studied earlier [^,|J by our 
methods. The results for the neutron cases, in particular, displayed K = ~ staggering not 
seen in the experimental results. We added Er to the mix because it has several K = \ bands 
and has been used in the past as a classic example of the Coriolis attenuation phenomenon. 

The special technical problems involved in applying our formalism have been fully doc- 
umented and solved in our previous papers and will not be repeated here. As far as the 
size of the model space is concerned we remind the reader that for the low energy data that 
we are analyzing, it suffices to describe the core by its ground-state band and those few 
excited bands, such as beta and gamma bands (but not restricted to these) that show some 
quadrupole collectivity with respect to the ground band. We then choose the size of the 
single-particle space coupled to these bands large enough so that further enlargement does 
not modify the results. In this respect we can claim to have essentially exact solutions of 
our semi-microscopic model, though the question of whether these are close to fully self- 
consistent solutions of the Kerman-Klein equations remains an open question. In general 
we have found that three major shells centered at the valence shell suffices. 

For all intents and purposes, our previous calculations contained only one free parameter, 
the strength n 2 of the quadrupole-quadrupole interaction, all other parameters having been 
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determined by other experimental results. (A minor cavil to this assertion is that small ad- 
justments of the spherical single-particle energies emerging from Woods-Saxon calculations 
were permitted for levels near the Fermi surface.) In the present work, we have a second 
parameter, the strength K\ of the magnetic interaction, to adjust. 

Before displaying the results found using the full theory, we illustrate the method of 
fitting the data with the additional parameter by studying a simplified model involving the 
coupling of a single j level to the core bands of Gd and with a value of the quadrupole 
interaction chosen so that we have a lower K — | band clearly separated from a higher 
K — | band. We see from Fig. [I| that for fixed k 2 , the staggering is reduced and finally 
eliminated as we raise K\ from its initial value of 0. At the same time the effect on the shape 
of the K = I band is negligible. Next, fixing the magnetic interaction at its optimum value, 
we see from Fig. |2] that changing the value of k 2 shifts band-head energies. What is not clear 
from these figures is that in the general case, the magnetic interaction also shifts relative 
band- head energies. Therefore in actual calculations, we first adjust k 2 to give the best 
results for other than the staggering, next find the optimum value of K\, and subsequently 
refit band heads with a new search on n 2 . This is an iterative procedure to find the best 
pair of coupling strengths. 

We next present the results of our calculations. The odd neutron nuclei 157 Gd and 159 Dy 
were studied previously in ||. We found that we could fit the energies satisfactorily except 
that a staggering of a K = 1/2 band for each nucleus that was not observed experimentally. 
To obtain the quality of agreement found in the previous work, it was necessary to include 
the P band, the 7 band, and an additional K = excited band. In Figs. |3] and f| we show 
the results of the augmented calculations and compare them to those found previously. It 
can be seen that the addition of the magnetic dipole interaction removes the unphysical 
staggering found earlier without any serious damage to the rest of the fit. 

The odd proton nucleus 157 Tb was used in a previous study of the status of the traditional 
core-particle model as an approximation to the KKDF model ||. In Fig. ^ we see that in 
this case no exaggerated staggering was predicted by the theory without a magnetic dipole 
interaction. The addition of a magnetic dipole-dipole field with strength of the same order 
of magnitude as for the other two nuclei gives, nevertheless, a perceptibly improved overall 
fit to the data. 

Last we have applied our theory to the case of 163 Er. We have studied both the negative 
and positive parity states of this nucleus. The results are shown in Fig. |6|. In the right 
panel, which displays the results for the levels of odd parity, one sees that we have found 
a good fit to the data, which includes three K — | bands. Remark that the fit is best for 
those two of these bands for which the data is more extensive. In the left panel, which is 
devoted to the positive parity levels, there is only a single K = \ bandhead. The lowest 
band, K — |, the only one that seems to be well-established experimentally, has features 
that we have failed to capture thus far. 

V. J J INTERACTION 

As mentioned in the introduction we can achieve results similar to those reported above 
by replacing the magnetic dipole interaction used previously by a term proportional to the 
square of the angular momentum. In Fig. |^ we contrast the new calculation thus engendered 
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for 157 Gd (right panel) with that found in the previous section (reproduced in the left panel. 
We find that both interactions have a similar effect on the staggering, but that the effect 
on the relative bandhead energies is different. In general, fits of comparable quality are 
obtained. To achieve the desired staggering, the strength, k[ of the JJ interaction should 
be set at K\ = 0.026. We get the best fit if the strength of the quadrupole field k 2 is altered 
slightly from what we had before: k 2 = 0.392 MeV/fm 2 compared to k 2 = 0.394 MeV/fm 2 . 

We compare the two interactions to understand their similarities and differences. If 
only matrix elements diagonal in both the core and the particle are considered, it is a well- 
known property of vector operators that the two interactions are equivalent. In fact this 
equivalence holds very well for the core, since for K = bands, with / values differing 
by two, the magnetic dipole operator has no off-diagonal elements. It is true that we lose 
this simplification for the gamma band, but the physics we are discussing is dominated by 
the ground band. When we turn, however, to the single-particle space, off-diagonal matrix 
elements of the magnetic moment operator can play an important role, as in generating 
magnetic dipole transitions, and are thus the agent responsible for the slightly different 
results found for the two interactions considered. 

VI. DISCUSSION 

The aim of this paper has been to present a reasonably convincing case, both on the basis 
of a priori arguments and by illustrative calculations, that we have uncovered the origin of the 
phenomena referred to as attenuation of the Coriolis coupling. By basing our considerations 
on a theoretical framework of great generality, we were able to rule out a large subset of 
explanations offered previously, since we obtained evidence of the attenuation phenomenon 
even when the physics associated with these explanations was naturally incorporated in our 
calculations. The only possibility remaining to us was to think about adding some term to 
the pairing and quadrupole interactions that define the Hamiltonian used in our previous 
work. Because of the nature of the effect sought, the most obvious choice is a magnetic 
dipole interaction. There are two such interactions that can be fitted easily into our semi- 
microscopic theory, the scalar product of the magnetic moment operator with itself and 
the square of the angular momentum operator. Calculations done with the former on four 
different nuclei showed that the excessive staggering encountered in previous calculations, 
but absent experimentally, could be made to disappear with a dipole interaction strength 
that remains constant within 10% for the nuclei studied. For one of the nuclei, 157 Tb, where 
no excessive staggering was calculated previously, we found a slight improvement to the 
overall fit with the same interaction strength. 

For one case, that of 157 Gd, we compared the results found for the two types of interaction 
mentioned above and found quite similar, but not identical, results. This suggests that 
we have identified the multipolarity of the additional interaction necessary to explain the 
quenching effect, but not its detailed form. 
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TABLES 



nucleus 


91 \Mu] 


ki [MeV//4] 


k 2 [MeV/fm 2 ] 


157 Gd 


0.371 


0.045 


0.394 


159 Dy 


0.345 


0.043 


0.368 


157 Tb 


0.362 


0.052 


0.420 


163 & (_) 


0.353 


0.04 


0.366 


163 Er (+) 


0.353 


0.04 


0.396 



TABLE I. Parameters used in the calculations reported in this work. The second column lists 
the gyromagnetic ratios of the core nuclei, the third column the values of the magnetic dipole 
interaction strenght, and the last the values for the quadrupole interaction strenght. 
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FIGURES 
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FIG. 1. A model calculation to illustrate the effect of the magnetic dipole-dipole interaction 
on staggering. A single-particle space with a single j value is coupled to the ground state bands of 
the even neighbors of 157 Gd. The strength of the quadrupole force is fixed at K2 = 0.564MeV/fm 2 . 
The strength of the magnetic dipole force is varied to demonstrate the effect of this interaction. 
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FIG. 2. The purpose of this figure is to demonstrate that the strength of the quadrupole field 
does not influence the staggering. The model is the same as in Fig [j]. The value, k\, of the magnetic 
interaction is set at the value that minimizes the staggering of the K = 1/2 band. Then the effect 
of varying is studied. 
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FIG. 3. Negative parity energy levels of 157 Gd. The calculation includes single particle states 
for three major shells and excited bands of the core. On the left is the previous calculation, with 
the strength of the quadrupole force set at 0.380 MeV/fm 2 and k\ = 0.0 MeV/^ 2 . On the right is 
the improved calculation with k 2 = 0.394 MeV/fm 2 and k\ = 0.045 MeV/^ 2 . 
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FIG. 4. Negative parity energy levels of 159 Dy. The calculation includes single particle states 
for three major shells and excited bands of the core. On the left is the previous result with the 
strength of the quadrupole force set at 0.383 MeV/fm 2 and k\ = 0.0 MeV/^ 2 . On the right is the 
improved calculation with k 2 = 0.368 MeV/fm 2 and «i = 0.043 MeV//i 2 . 
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FIG. 5. Positive parity energy levels of 157 Tb. The calculation includes single particle states 
for three major shells and excited bands of the core. On the left is shown the previous result where 
the strength of the quadrupole force is 0.428 MeV/fm 2 and k\ = 0.0 MeV/^u 2 ,. On the right is the 
improved calculation with k 2 = 0.420 MeV/fm 2 and k\ = -0.052 MeV/V 2 ,. 
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FIG. 6. Energy levels of 63 Er. Single particle states for three major shells and excited bands 
of the core are included. On the left are the positive parity states and on the right the negative 
parity states. For the positive parity states the strength of the quadrupole force is 0.413 MeV/fm 2 
when ki = 0.0 MeV/^x 2 ( no t shown) and k 2 = 0.399 MeV/fm 2 when k\ = -0.040 MeV/^ 2 . For 
the negative parity states the strength of the quadrupole force is 0.375 MeV/fm 2 when k\ = 0.0 
MeV/V 2 (not shown) and k 2 = 0.366 MeV/fm 2 when ki = -0.040 MeV//i 2 . 
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FIG. 7. Comparison of two magnetic interactions. On the left is the same fit as in Fig. || using 

the magnetic dipole operator and on the right the fit using the angular momentum form. For the 

calculations in the left panel, we have k 2 = 0.394 MeV/fm 2 with k x = 0.045 MeV/// 2 ;. for the 
calculation on the right, k 2 = 0.392 MeV/fm 2 and k[ = 0.026 MeV//i 2 . 
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